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A fraction of the thermal protons in the outer envelope of an X-ray cluster have velocities that 
exceed the local escape speed from the cluster gravitational potential. The Coulomb mean-free- 
path of these protons is larger than the virial radius of the cluster at temperatures > 2.5keV. The 
resulting leakage of suprathermal particles generates a coUisionless shock in neighboring voids and 
fills them with heat and magnetic fields. The momentum flux of suprathermal particles cannot be 
confined by magnetic tension at the typical held strength in the periphery of cluster halos ^^G). 
Over a Hubble time, thermal evaporation could drain up to a tenth of the cluster gas at its virial 
temperature. The evaporated fraction could increase dramatically if additional heat is deposited into 
the gas by cluster mergers, active galactic nuclei or supernovae. Thermal evaporation is not included 
in existing cosmological simulations since they are based on the fluid approximation. Measurements 
of the baryon mass fraction in the outer envelopes of hot clusters (through their Sunyaev-Zel'dovich 
effect or X-ray emission) can be used to empirically constrain their evaporation rate. 



I. INTRODUCTION 

It is commonly assumed that the fluid approximation 
applies to the virialized gas in X-ray clusters. However, 
the Coulomb mean-free-path of a proton or an electron 
in a plasma of temperature T and density contrast A 
relative to the present-day (z = 0) mean cosmic density 
is given by, 



A(0) = 1.3Mpc(M;)' 
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A value of A ~ 25 characterizes the local gas density at 



the virial radius 



of cluster dark-matter halos 



For a cluster that virializes at the present time 
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where we adopt the standard set of cosmological parame- 
ters 'i'l and assume primordial composition with 24% he- 
lium by mass. Thus, the mean- free-path of protons and 
electrons is comparable to the size of the virialized en- 
velope of hot X-ray clusters. Hydrodynamic simulations 
1^130) 13 show that a substantial portion of the surface of 
these envelopes is in contact with rarefied voids in which 
A < 1. Protons moving away from the cluster envelope 
towards a neighboring void must be confined by gravity 
since Coulomb collisions in the surrounding A < 1 envi- 
ronment are too rare to confine them. This applies even 
better to protons with a kinetic energy E larger than the 
mean thermal value, since the Coulomb mean-free-path 
scales as E"^ . If the kinetic energy of a proton exceeds the 
local gravitational binding energy, gravity cannot keep 
the proton bound to the cluster. Thermal plasma would 
therefore evaporate from the cluster (as each escaping 
proton must be accompanied by an electron so as to pre- 
serve charge quasi- neutrality) . This process reflects a de- 
parture from the fluid description of the intracluster gas 



and cannot be captured by numerical hydrodynamic sim- 
ulations which are based on the fluid approximation. The 
purpose of the this paper is to gauge the significance of 
this evaporation process in X-ray clusters and to motivate 
future numerical simulations involving kinetic theory in 
cosmological three-dimensional geometries. Significant 
evaporation would modify the baryon fraction in clus- 
ters and weaken the case for it having a universal value 
that can be used for cosmological distance determination 
0, 1^ Hol Ull | 3 p . The evaporation would also heat the 
gas in voids 13* and potentially change the Lya absorp- 
tion signature of voids in quasar spectra at low redshifts 
. The fractional helium abundance in clusters should 
also increase as only hydrogen evaporates. 

The use of X-ray clusters for cosmological parameter 
determination (e.g. regarding the dark energy equation 
of state) is one of the active frontiers in cosmology. A 
large number of surveys for the Sunyaev-Zeldovich effect 
as well as the X-ray emission of clusters are currently 
being conducted. The observational data is usually com- 
pared to cosmological hydrodynamical simulations. In 
this paper we point out that the fluid approximation is 
not valid in the outer parts of hot clusters, where thermal 
evaporation of gas is likely to take place. The evapora- 
tion needs to be included in future studies attempting to 
refine the use of X-ray clusters for precision cosmology. 



II. MASS LOSS RATE 

The escape velocity Wcsc from a virialized cluster is de- 
fined as the velocity necessary to overcome the gravita- 
tional potential energy per unit mass W of the cluster, 
5^csc = —W{y). The virial theorem relates the mass- 
averaged value of W to twice the mean kinetic energy 
per unit mass in the cluster {—W) = {v"^) and implies 



(^csc) ~ 4(t!^), where angular brackets denote mass- 
averaged values . 

In the envelopes of clusters, the timescale for equili- 



2 



bration of the electron and proton temperatures through 
Coulomb collisions is comparable to the Hubble time 
|l6Lll7l |. However, within any plasma confined by a grav- 
itational field, there is an electric field which balances the 
pressure gradient of the electrons |0| , 



E 



(3) 



where Pe = UeksTe is the electron pressure and e is the 
electron charge. This electric field binds the electrons, 
which are otherwise too light to be bound by gravity. The 
same field reduces the depth of the gravitational barrier 
for a proton to escape. In an electron-proton plasma with 
a single temperature, the gravitational potential barrier 
is reduced by a factor of 2. If Te = then E — and the 
barrier equals the full gravitational potential. The evapo- 
ration rate depends only on the mean thermal energy per 
electron-proton pair irrespective of how they share this 
energy. For simplicity, we will assume that all particles 
have the same temperature from now on. 

The local barrier in the outer envelope of the cluster 
is smaller than its mass-averaged value across the en- 
tire cluster. The evaporation rate depends on the local 
value of the parameter, 77(r) = v'^^^/v'^, at the envelope. 
For a Maxwellian velocity distribution of an equilibrium 
temperature T, the average kinetic energy per proton is 
impf^ = ^ksT, where rrip is the proton mass. The frac- 
tion e(r) of the thermal particles that exceed the escape 
speed is given by 



3>)/2 



exp(-C)C'/'rfC- 



(4) 



For 1 < r; < 4, we find e « lAle''^-^'^ to better than 10%. 
As particles from the tail of the Maxwellian distribu- 
tion escape the cluster, this tail is re-populated through 
Coulomb collisions in the higher density gas bordering 
the envelope from within the cluster. 

Next we calculate the fraction of protons above the es- 
cape speed as a function of fractional radius R — (r/rvii ). 
Hydrostatic equilibrium for a primordial composition im- 
plies 



kBT{r) 



1.18 



GM{r')p^{r') 



dr', 
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and the escape speed for a proton (after subtracting off 
the repulsive electric force from the gravitational force) 



= M = 0.82 



GM(r') 



dr', 



(6) 



where Pg{r) is the gas mass density and M{r) is the total 
interior mass at radius r. Numerical simulations and 
X-ray data indicate that the hot gas follows the dark 
matter distribution in the outer parts of X-ray clusters 



[l^ . We therefore adopt the NFW 0| density profile, for 
which M{x) oc [ln(l -I- a;) - x/{l + x)] and pg cx l/[a;(l + 
x^], where x ~ r/rg and rg = r^ij-fC with C being the 
concentration parameter. The normalization of the virial 
mass or virial radius scale out from the expression for 
1] — v1^^/[?>kBT /rrip], which depends only on R and C. 
The radial profile of the gas temperature is calculated 
from Eq. ©. 

Figure 1 depicts e(i?) for three values of C. We use 
the integration upper limit of r' = lOrvir instead of oo in 
the above integrals, as it represents the typical radius of 
a void or half the distance between rich clusters. For the 
characteristic value of recently collapsed halos, C ^ i, 
the esca pe f raction is ~ 7.6% at i? ~ 1, corresponding to 
rj ^ 2.3 |42|. The asymptotic power-law slope of —3 for 
the radial density profile yields 77 — 1.87 and e « 13% for 
large values of C. Since pg oc R~^ at R ^ 1/C, fresh 
suprathermal particles can be supplied to the coUisionless 
envelope from the highly coUisional interior at i? < 0.5. 

Clearly, the cluster plasma has to be sufficiently hot in 
order for the suprathermal particles not to be reflected 
back into the cluster through Coulomb collisions. Their 
mean- free-path (oc rj^) exceeds the virial radius of a clus- 
ter for temperatures 



k, 



.T>2.5keVx(l + .)(^) - 
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where we used Eqs. ^ and Q with the additional scal- 
ing to arbitrary redshift, z. 

Hydrodynamic simulations 0, 0, 01 indicate that 
infalling gas is often channeled into clusters through fila- 
ments, while other parts of the cluster border low-density 
voids. Thermal evaporation will occur preferentially 
across the contact surface with voids where there is negli- 
gible resistance to the outward flux of escaping particles. 
Magnetic fields which are tangled on small scales would 
suppress evaporation since the proton Larmor radius is 
negligibly small. Such fields could ori gina te either from 
galactic winds 20] and quasar outflows |2lLl22l| or be gen- 
erated in the coUisionless shocks that virialize the cluster 
gas ji^ 1^ ■ The fields are likely to be accompanied by 
relativistic cosmic-rays which are produced by the same 
sources. Indeed merging X-ray clusters are known to be 
surrounded by radio halos due to synchrotron emission 
by relativistic electrons [23ll25| . The cosmic-ray pressure 
mig ht subject the magnetic field to the Parker instability 
|2(tI |. producing channels of open field lines along which 
thermal evaporation will not be inhibited. In deriving the 
evaporation rate below, we include a suppression factor 
/ < 1 to account for the fraction of the cluster surface 
area where evaporation is not suppressed by inflowing gas 
or by magnetic fields. 

The mass evaporation rate of cluster plasmas that sat- 
isfy Eq. Q can be estimated from the outward flux of all 
thermal particles above the escape speed of the cluster 



3 



0.1 



O.OB - 




0.2 0.4 0.6 0.8 



FIG. 1: Fraction e of protons above the escape speed as a func- 
tion of fractional radius R = (r/rvir) for an NFW Q density 
profile with different values of the concentration parameter 
C = 2 (dashed line), 4 (dotted) and 8 (solid). 



at the virial radius l5C 



Mg = 47rx0.76/epg(rvir)r^irWosc ~ Sferj 



1/2 
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where Mg is the interior gas mass at the virial radius and 
t = H{z)~^ is the Hubble time at the redshift of interest. 
Hence, the fraction of the cluster gas that evaporates over 
a fraction a of the Hubble time is ~ Sa/e??^/^. For af ~ 
1/3 and rj ~ 2.3, the total evaporated fraction is ^ 10%, 
comparable to the fraction of baryons which are locked in 
stars [27!]. A deficit larger by a factor of a few (attainable 
with / ~ 1 and 77 ~ 2) is still allowed by existing data 
on the baryon mass fraction in cluster envelopes 12, 0, 
im . As accretion will taper off in the future of a A- 
dominated universe 29, "3^, the hot envelopes of X-ray 
clusters might be depleted significantly over a timescale 
> 10"yr- 

In order for magnetic confinement to operate, the mag- 
netic tension B^/Stt needs to exceed the momentum fiux 
of the suprathermal particles ~ ePg^'escI ^Ise the mag- 
netic field will be moved around by the suprathermal 
wind. Although the field is anchored to the bulk of the 
plasma, a large load of suprathermal particles will re- 
shape its topology and allow escape, e.g. through the 
"garden-hose" instability [U I13and other buckling in- 
stabilities of the magnetic field [23 that will open chan- 
nels of outflow for the suprathermal particles. These in- 
stabilities are sourced by the anisotropic distribution of 



particles in the local plasma frame and are particularly 
powerful when the pressure anisotropy of the plasma ex- 
ceeds the magnetic pressure. The required field strength 
to resist the outward momentum fiux of suprathermal 
particles is 



B > ImG 



2.3 



5 kcV 



1/2 



(9) 



Synchrotron and hard X-ray observations as well as 
Faraday rotation measurements in the periphery of the 
extended halos around clusters infer a magnetic field 
strength which is 3-10 times lower than the required mag- 
nitude 0, . The associated magnetic tension is 
at least an order of magnitude smaller than required and 
could therefore allow significant leakage of suprathermal 
particles. 

The mean-free-path of suprathermal protons might be 
reduced through scattering on plasma waves [371 Issf . 
which are excited in the intracluster plasma by the col- 
lisionless shock formed when accreting gas impacts on 
the cluster boundary or through the turbulence gener- 
ated by the motion of sub-halos within the cluster. How- 
ever, the excitation of plasma turbulence is expected 
to be unsteady and the waves would leak out during 
quiescent episodes of weak accretion. The long mean- 
free-path of protons makes a large viscosity coefficient 
that would have suppressed acoustic turbulence alto- 
gether in the absence of magnetic fields. (Although ex- 
isting hydrodynamic simulations identify acoustic tur- 
bulence in cluster cores, their limited numerical resolu- 
tion and incomplete treatment of plasma physics do not 
allow them to examine reliably the turbulence damp- 
ing rate Plausible excitation mechanisms of 
Alfven waves in a relaxed cluster produce negligible 
power on wavelengths that resonate with the microscopic 
ion gyro-radius, ~ 10^^ cm(_B/0.1/iG)^^; under such 
circumstances the waves cannot scatter effectively the 
suprathermal particles and suppress their large momen- 
tum flux. 

The outgoing beam of suprathermal particles would 
generate a shock in surrounding voids through the plasma 
two-stream instability j3ll |. The dispersion relation 
of the instability at finite temperatures implies that 
a Maxwellian-tail beam will not excite this instability 
within the hot cluster plasma but only as it enters the 
much colder environment of a void. The resulting shock 
will heat the gas in the voids and enrich it with magnetic 
fields. 



III. DISCUSSION 

We have shown that if magnetic confinement is inef- 
fective (Eq. O, then thermal evaporation is capable of 
removing up to a tenth of the cluster gas at its virial 
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temperature over a Hubble time (E g. |Ht. Local en- 
ergy injection from cluster mer gers |41| . active galac- 
tic nuclei or supernovae l43| could enhance the 
evaporated fraction as it would raise the gas tempera- 
ture and lower the corresponding value of rj. The evapo- 
rated fraction is exponentially sensitive to enhancements 
in the gas temperature, with e (x exp(— 1.377). For a fixed 
energy injection per baryon (as expected from a con- 
stant star/black-hole formation efficiency), the change in 
7] would be largest at low values of T. Indeed, the baryon 
fraction in the shallower potential wells of poo r clusters 
appears to be smaller than in rich clusters j23|, indicat- 
ing gas loss through additional heat supply to the virial 
energy budget. 

An essential ingredient for the evaporation to be main- 
tained is that the tail of suprathermal protons will be 
replenished despite the low coUisionality of the plasma 
in the cluster atmosphere. Indeed protons just below 
the escape speed are turned around by gravity and then 
fall back towards the cluster center. Along their trajec- 
tory they encounter higher density gas, where they collide 
with other protons and re-populate the Maxwellian tail 
once again. An approach to the coUisional regime merely 
requires sinking by a factor of 2 in radius for the NFW 
density profile, where pg cx for r^h/C ^ r < rvw- 
This replenishment process does not occur for weakly in- 
teracting dark matter particles, which would therefore 
not evaporate steadily in relaxed clusters. The dark mat- 
ter might nevertheless experience transient evaporation 
episodes as a result of time-dependent gravitational in- 
teractions during violent mergers. 

Although existing simulations of clusters are performed 
in the fluid approximation, they can nevertheless be im- 
plemented for calculating the thermal evaporation rate 
under more realistic conditions than those considered in 
Fig. 1. The simulated distributions of the gravitational 
potential and the gas temperature can be used to derive 
Wosc(r) and e(r) within the region where the plasma is 
coUisionless. A corresponding evaporation term can then 
be added to the local conservation laws. The significance 
of magnetic confinement could be examined through 
magneto- hydrodynamic simulations with the im- 

portant inclusion of the stress tensor for suprathermal 
particles and cosmic-rays. 

Better observational constraints on the intracluster 
magnetic field [3 HE HE QEl could determine whether 
the inferred magnetic tension in cluster halos is indeed 
unable to sustain the momentum flux of suprathermal 
protons above the escape speed (Eq. EJ. The same 
fields affect intracluster cosmic-ra ys w hich produce neu- 
tral pions that decay into 7-rays [43 ■ Future 7-ray ob- 
servations with GLAST^il'l could constrain the magnetic 
confinement time of cosmic-rays, and by extrapolation - 
suprathermal particles. 

Exis ting X-ray data on the baryon mass fraction 
0, Hol l28| probe only the interior of clusters {R < 0.3), 



where A < r. The observed baryon fraction there is 
already known to be smaller than the expected cosmic 
average. Future extensions of existing data sets to the 
low-density envelopes of hot clusters through deep ob- 
servations of their X-ray emission or Sunyaev-Zel'dovich 
effect could determine whether there is a deficit in 
the baryon fraction relative to the cosmic value beyond 
the level expected from the consumption of gas by star 
formation 
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